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The transfer of weakly bound KRb molecules from levels just below the dissociation threshold into 
the vibrational ground state with shaped laser pulses is studied. Optimal control theory is employed 
to calculate the pulses. The complexity of modelling the molecular structure is successively increased 
in order to study the effects of the long-range behavior of the excited state potential, resonant spin- 
orbit coupling and singlet-triplet mixing. 
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I. INTRODUCTION 

Research on cold and ultracold molecules has been one 
of the most active areas of atomic and molecular physics 
over the last decade and continues to draw much atten- 
tion [lj . Current activities are by and large still focussed 
on producing (ultra)cold molecules, either by direct cool- 
ing or by assembling molecules from cooled atoms using 
external fields. While direct methods such as Stark de- 
celeration [2| have not yet reached the regime of ultra- 
cold temperatures (T < 100 /xK), photo- and magneto- 
association create molecules in their electronic ground 
state that are ultracold but in highly excited vibrational 
levels [3-[a|. However, prospective applications ranging 
from high-precision measurements to quantum informa- 
tion carriers [l| require stable ultracold molecules. This 
has triggered the quest for molecules in their absolute 
ground state. A major breakthrough toward this long- 
standing goal was achieved when several groups produced 
molecules in the lowest rovibrational level of an electronic 
ground statepotential via Stimulated Raman Adiabatic 
passage @, ml, p hotoassociation [8|, Q and vibrational 
laser cooling [lfj. Finally, the ability to control not only 
the rovibrational but also the hyperfine degree of frccdonr 
has recently paved the way toward ultracold molecules in 
their absolute ground state [Tlj . 

An earlier, alternative proposal to reach molecules 
in their vibronic ground state was based on employing 
laser pulse shaping capabilities: An optimally shaped 
laser pulse can coherently transfer, via many Raman 
transitions, vibrationally highly excited molecules into 
v = 12]. The attractiveness of this proposal rests on 
the fact that optimization is carried out iteratively, both 
in experiment and in theory, and does not require de- 
tailed knowledge of the molecular structure to identify 
optimal pulses. One might nevertheless ask whether and 
how qualitative changes in the molecular structure affect 
the optimal solution. This question defines the scope of 
the present work. We solve the same optimization prob- 
lem - population transfer from a weakly bound vibra- 
tional level just below the dissociation threshold of the 
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electronic ground state potential to v = - for two differ- 
ent molecules, KRb and Na2, successively taking effects 
into account that strongly alter the molecular structure 
of KRb. Our choice to focus on this molecule is motivated 
by the long-standing dEMl] and continuing [H QjJ HJJ 
experimental efforts on KRb. 

The paper is organized as follows: The general theoret- 
ical approach is described in Section [TTJ Sections IIIUlVl 
present and analyze the short shaped laser pulses that 
achieve the vibrational transfer. The complexity of the 
model for the molecular structure of KRb is successively 
increased, cf. Fig. [TJ Based on a two-state model, Sec- 
tion |nT] compares the optimization of vibrational transfer 
for KRb and Na2 molecules. They differ in the long- 
range behavior of their excited state potentials, 1/R 3 
vs 1/R 6 . Section HVl is dedicated to the effect of spin- 
orbit interaction in the electronically excited state which 
may lead to resonant coupling between a singlet and a 
triplet state |2ll - l23J . This is investigated with a three- 
state model in Section IIV1 Adding one more channel, 
Section [V] studies the transfer of a coherent superposi- 
tion of singlet and triplet molecules in their electronic 
ground state into the ground vibrational level of the sin- 
glet potential. This is possible for polar molecules due 
to the broken gerade-ungerade symmetry Q . Section IVII 
concludes. 



II. THEORETICAL APPROACH 

A. Model 

The linear Schrodinger equation describing the in- 
ternuclear dynamics of two atoms is considered. For 
molecules formed in a quantum degenerate gas, the 
many-body dynamics are then neglected. This approach 
is justified by the time scales of standard optical and/or 
magnetic traps 24] . While the internuclear dynamics and 
pulse shaping occur on the time scale of picoseconds, the 
many-body dynamics for conventional traps is character- 
ized by microseconds. The many-body system will have 
to adjust to the new internal state, but this is going to 
happen on a much slower time scale only after the pulse 
is over 123. 
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FIG. 1: (color online) Potential energy curves and initial and target wave functions for vibrational stabilization of Na2 (left) 
and KRb (right) molecules. The complexity of the model for the molecular structure of KRb is successively increased from two 
to four channels. 



The time-dependent Schrodinger equation, 

»^Wt)) = (Ho + /te(t))Mt)>, (1) 

is solved with a Chebychev propagator [26[. The Hamil- 
tonian comprises two or more electronic states as spec- 
ified below, cf. Fig. [TJ The interaction with the laser 
field is treated in the dipole approximation. The rotating 
wave approximation is not invoked in order to allow for 
strong fields and multi-photon transitions. The Hamil- 
tonian and the wavefunctions, tpi(R) — (R\tpi), where i 
denotes the channel, are re pres ented on a Fourier grid 
with adaptive grid step [27H29|. In order to obtain vi- 
brational eigenfunctions and Franck-Condon factors, the 
Hamiltonian, Hq, is diagonalized. 



B. Optimal control theory 



to use dynamical information from intermediate times, a 
new functional is defined, 



J 



-F- 



g(e)dt. 



(4) 



where g(e) denotes an additional constraint over the field. 
Often g(e) is chosen to minimize the pulse fluence. This 
implies a replacement rule in the control equation for the 
field. However, a choice of g(e) that leads to vanishing 
changes in the field as the optimum is reached may be 
advantageous from the point of view of convergence [30| . 
It is employed here, 
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Denoting the formal solution of the Schrodinger equa- 
tion at time t by 



\i>(t)) = u{t, o)|V(o)> , 



(2) 



the objective functional for a transition from initial state 
\ipini) to target state \iptarget) at the final time T can be 
written 



F= \(^ini\yi (T,0;e)\i; target }\ 



(3) 



It corresponds to the overlap of the initial state that has 
been propagated to time T under the action of the field 
e{t) with the target state. A field is optimal if it com- 
pletely transfers the initial state, \ipmi}, to the target 
state, \tptarget), i-e. if F reaches a value close to one. 

The objective F is a functional of the field e(t). It 
explicitly depends only on the final time T. In order 



and e(t) is taken to be the field of the previous iteration. 
This choice of g(e) implies an update rule rather than 
a replacement rule in the control equation for the field. 
Physically, it corresponds to minimizing the change in 
pulse energy at each iteration. The shape function S(t), 
S(t) = sin 2 (nt/T), enforces a smooth switch-on and off 
of the field. The parameter A controls the optimization 
strategy: A small value results in a small weight of the 
constraint, Eq. ([5]), allowing for large modifications in the 
field, while a large value of A represents a conservative 
search strategy with only small modifications in the field 
at each iteration. 



We employ the Krotov algorithm [30-34] to obtain the 
control equations. The derivation of the algorithm for the 
target, Eq. ([3]), and the constraint, Eq. (JSJ) , is described 
in detail in Ref. [l2(. It yields the following prescription 
to improve the field, 



e 3 (t) 



i(t) 



s{t) n 



Ulpini\V {T,0; £j-l)\tptarget){^target\(i (*, T; £j_i) p, 0(t, 0; £,-)|V>mm) } 



(6) 



r 



at the jth iteration step. The first overlap in the paren- 
thesis can be shown to be time- independent. The second 
overlap contains a backward propagation of the target 
state from time T to time t under the old field, £j_i(£), 
and a forward propagation of the initial state under the 
new field, Sj(t). Eq. © is implicit in £j(t). This is reme- 
died by employing two different grids in the time dis- 
cretization, see Rcf. f 2] for details. 



III. THE ROLE OF THE LONG-RANGE 

BEHAVIOR OF THE EXCITED STATE 
POTENTIAL: COMPARING KRB TO NA 2 

The investigation is started with a simple two-state 
model. The Hamiltonian reads 



c ,, f + l W,( R ) 



£e(i) 



T + V A1 



A l S7 



(R) 



(7) 



where (g/u) applies only to Na2. T denotes the kinetic 
energy operator and V^(R) the potential energy curve 
of channel i. For the Na2 molecule, we employ the 
same potential energy curves as in Ref. Il2lll54l which 
were obtained from molecular spectroscopy |35l . |36| . The 
ground state potential for KRb is also known spectro- 
scopically [37| . The excited state potential energy curves 
of KRb are unfortunately not yet known with spectro- 
scopic precision. We have therefore employed the re- 
sults of ab initio calculations [38| at short internuclear 
distances together with an asymptotic expansion of the 
form V asy (R) = — -j& — ^f with the long-range coefficients 
taken from Ref. [39[. The laser field, e(t), couples to the 
molecule via the transition dipole moment, fi. The lat- 
ter is expected to depend on internuclear distance at least 
for short R. However, it is approximated by /t = I a.u. 
since no data on the -R-dependence could be found in the 
literature. 

The initial state, ipini(R) = (R\4>ini) — ^P v g {R) is taken 
to be a vibrational wavefunction of the electronic ground 
state corresponding to a weakly bound level just below 
threshold. Counting the levels downward from the last 
bound one, v = vi as t — 3 for Na2 and v = vi ast — 6 for 
KRb. These levels were chosen to yield a comparable 
binding energy for the two species. They are very simi- 
lar to levels that would be populated when molecules are 
formed from atoms using a Feshbach resonance or pho- 
toassociation but with a somewhat larger binding energy. 
This approximation eases calculations since, as discussed 
below, the absolute value of the optimization time is de- 
termined by the binding energy of the initial state. It 



does not affect the comparison between Na2 and KRb 
which is based on relative times. 

While the model, Eq. ([7]), does not capture the full 
complexity of the molecular structure, it allows for a 
straightforward comparison with earlier work on the Na2 
molecule [12|]. In particular, it is useful to highlight the 
influence of the long-range behavior of the excited state 
potential which is 1/i? 6 for KRb but 1/R 3, for Na2- From 
a spectroscopic point of view, one expects the 1/R 6 - 
behavior of KRb to be more favorable for vibronic tran- 
sitions between highly excited vibrational levels. This 
is attributed to better Franck-Condon overlaps between 
alike potentials. However, from a dynamical point of 
view, a maximum difference potential between electronic 
ground and excited state, AV(R) = ^(V e (R) — V g (R)), is 
more desirable. A larger difference potential is obtained 
for an excited state with 1/R 3 long-range behavior, i.e. 
for Na2- The reasoning behind this argument is that the 
wave packet that launched by the pulse on the excited 
state experiences a much larger gradient and can better 
accelerate its motion toward shorter distances |40l |. 

In order to decide which of the two arguments is rele- 
vant for vibrational stabilization of ultracold molecules, 
we determine the minimum pulse energy and minimum 
optimization time required to achieve a transfer to the 
target state of 99% or better. This is based on the fact 
that many solutions to the control problem exist. Which 
solution will be found by the algorithm depends cru- 
cially on the boundary conditions - optimization time 
and pulse power or pulse energy. However, if the re- 
sources in terms of time and energy are insufficient, no 
solution will be found. The lower limits to T and £p can 
thus be used to characterize the control problem and the 
solution strategy. 

Figure [2] compares the optimal pulses and their spec- 
tra for Na2 (left) and KRb (right). These pulses are the 
results of a three-step optimization. Initially, an opti- 
mization time of T = 16 ps corresponding to about twice 
the vibrational period of the initial state of Na2 was cho- 
sen [12J. The guess pulses for this optimization were 
constructed as a series of 100 fs pulses with two central 
frequencies reflecting the peaks of the Franck-Condon 
factors of the initial and the target states. This choice 
provides a large enough spectral bandwidth for given T. 
The pulses were optimized until a transfer of 99% or bet- 
ter was achieved. For the second optimization step, the 
pulses were compressed in time following the recipe de- 
tailed in Ref. [liq . i.e. points e(wj) were removed such 
that Auj is increased and T decreased. The resulting 
pulses were then employed as guess pulses for the sec- 
ond step of optimization. For Na2, compression of T by 
a factor of 8 was possible while for KRb a factor of 4 
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FIG. 2: Optimal pulses (upper panel) and their spectra (lower 
panel) for Na 2 (left) and KRb (right). 



turned out to be the limit. For larger compression fac- 
tors, the optimization did not result in any appreciable 
population transfer to the target state. The difference 
in the minimum optimization time T for Na2 and KRb 
is explained in terms of the different reduced masses of 
the molecules. Although the binding energy of the initial 
states is comparable, cf. Fig.[TJ the motion of the heavier 
KRb is slower. This is reflected in the vibrational period 
of the initial states, 8ps for Na2 compared to 14 ps for 
KRb. Therefore, the compression factor taken with re- 
spect to the corresponding vibrational period, is 4 for 
Na2 compared to 3.5 for KRb, i.e. the maximum fac- 
tor for compression in time is very similar. Finally, the 
minimum energy required for optimal transfer was deter- 
mined in step three where the guess pulses were taken 
to be the optimal pulse of the previous step divided in 
amplitude by some factor. If the factor was too large, op- 
timization did not result in any appreciable population 
transfer, otherwise 99% transfer or better were achieved. 
This way, a sharp limit for the mininum required pulse 
energy was obtained. It amounts to 78 /iJ for Na2 and 
61 /./J for KRb where the focal radius of the laser beam 
is assumed to be 100 tim. 

The lower bounds on the optimization time and the 
pulse energy represent the main difference between the 
optimal pulses for Na2 and KRb. As can be seen in Fig.[2j 
the overall temporal and spectral structures of the pulses 
are fairly similar. The spectral bandwidth of the opti- 
mal pulse is significantly larger for Na2 than for KRb. It 
is, however, difficult to attach a physical meaning to the 
spectral bandwidth. We do not filter out undesired spec- 
tral components [4l|,[42|], i.e. those components that do 
not correspond to given vibronic transitions in order to 
avoid significantly increased convergence times and nu- 
merical effort. A constraint formulated directly in fre- 
quency domain to contain the bandwidth within a cer- 
tain spectral reg ion cannot be enforced simultaneously 
with Eq. ||2J) 12]. Therefore, with progressing iteration, 
the spectral bandwidth of the optimized pulses grows. 
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FIG. 3: (color online) Projection of the time-dependent wave 
packets onto the vibrational eigenstates of the electronic 
ground (bottom) and excited (top) states for Na 2 (left) and 
KRb (right) 
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FIG. 4: (color online) Franck-Condon factors of the ground 
state levels with all excited state levels for Na 2 (left) and KRb 
(right). The dissociation limit of the electronic ground state 
defines the zero of energy and the dissociation limit of the 
excited state is 16965 cm -1 for Na 2 and 12737 cm -1 for KRb. 



Since in principle this growth can be suppressed by fil- 
tering [4l|, [42| , it is rather an artifact of the algorithm 
than a physically significant finding. 

The dynamics induced by the optimal pulses are anal- 
ysed in Fig. [3] by projecting the ground and excited state 
components of the time-dependent wave packets onto 
the vibrational eigenstates. The dynamics show simi- 
lar features for Na2 and KRb. The optimal pulse de- 
pletes the initial state and pumps most of the population 
to the electronically excited state, distributing it over a 
wide range of highly excited vibrational levels. The wave 
packet then climbs down the potentials. In the middle 
of the optimization time interval, it reaches ground state 
levels v w 30 (v ~ 40) for Na2 (KRb) ground state which 



is about half way down to the bottom of the potential 
well. In the last quarter of the optimization time inter- 
val, population is accumulated in a superposition of a few 
excited state levels around v' fa 10 for Na2 (v' fa 18 for 
KRb). In a final step, this superposition is transfered to 
the target level, v = 0. 

The population dynamics is rationalized by an analy- 
sis of the Franck-Condon map, cf. Fig. @] The dynamics 
start in the upper right corner of the map and follow 
the main ridge until v « 30 (v w 40) for Na 2 (KRb) 
is reached. At this point the Franck-Condon map takes 
approximately the shape of a parabola where the right 
branch connects to v fa 30 (v fa 40) Na 2 (KRb) while the 
left branch connects to the target level, v = 0. This ex- 
plains the dynamics in the second half of the optimization 
time interval where the population is pumped into the ex- 
cited state levels, v' ~ 10 for Na 2 (v' fa 18 for KRb), that 
are reached by the two branches of the parabola, i.e. that 
are the ideal gateway to v = 0. 

To conclude the comparison of Na 2 and KRb, the min- 



imum optimization time is dictated by the mass of the 
molecule and population transfer in terms of required 
pulse energy is more favorable for a 1/R 6 than a 1/P 3 
excited state potential. The overall dynamics of the pop- 
ulation transfer is rather similar for the two molecules 
and is easily rationalized by the structure of the Franck- 
Condon map. 



IV. SPIN-ORBIT COUPLING IN THE 
ELECTRONICALLY EXCITED STATE 



The complexity of our model is increased to take spin- 
orbit interaction in the electronically excited state of 
KRb into account. Similarly to the Rb 2 and RbCs 
molecules, spin-orbit interaction may lead to resonant 
coupling and strong non-adiabatic effects |21H23l |. This 
is captured by a three-state Hamiltonian, 



J 



A 3 ,(t) = 



/ T + y xls+ (R) £e(i) 

Ae*(t) t + VW(R) W™(R) 

Wf£(R) t + lW(R)-Wj£(R) y 



(8) 



r 



where the diagonal and off-diagonal spin-orbit interaction 
terms, W 7 ™ (R) and W™{R), are introduced. In prin- 
ciple, the 6 3 II excited state has a dipole coupling with 
the lowest triplet state. One can thus transfer molecules 
from the lowest triplet state via Raman transitions into 
the singlet ground state. This will be investigated be- 
low in Section [Vj and the present section is devoted to 
studying the effect of non-adiabaticities in the electron- 
ically excited state. Here, the initial state is purely sin- 
glet, i.e. the same weakly vibrational level of the X l Y, + 
state, v = 93 = vi as t — 6, as in the previous section is 
considered. 

Unlike in the case of Rb 2 where the spin-orbit in- 
teraction terms were determined spectroscopically [43j . 
no such accurate data is available for KRb. We have 
therefore resorted to the parametrization of W™ (R) and 
WgQ (R) in terms of Morse functions, 
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that was introduced by Bergeman et al. for RbCs j44j . 
These two functions show a dip at intermediate distances 
and level off toward a constant value at long range. As 
a first guess, we have emp loyed the values for the pa- 
rameters P/ from Ref. [44j scaled to reproduce the cor- 
rect asymptotic limit of the fine-structure splitting of 
237.595 cm -1 . The corresponding values for the param- 



eters Pi are listed in the first row of Table |U Since the 
parameters P/ are not accurately known, we have varied 
the Pi in order to estimate the maximum effect that the 
spin-orbit interaction can have on the vibrational wave 
functions and Franck-Condon matrix elements. This pro- 
vides the starting point for studying the strongest possi- 
ble effect of the spin-orbit interaction on the optimization 
and the dynamics under the optimal pulse. Two differ- 
ent choices of spin-orbit coupling are employed, referred 
to below as cases 1 and 2. The corresponding parame- 
ters are listed in the second and third rows of Table |U 
The modification of the parameters is quite substantial 
and larger than what can realistically be expected. How- 
ever, the point here is to demonstrate the most positive 
and most negative effect that the spin-orbit coupling may 
have on the vibrational stabilization and to explore its 
influence on the optimization. 

In case 1, we have modified Pp n and P™. The latter 
corresponds to the width of the dip in the off-diagonal 
spin-orbit coupling, while the former represents the con- 
stant offset of the diagonal spin-orbit coupling, which es- 
sentially causes a relative shift of the vibrational ladders 
of the A^T, and 6 3 n states. This choice of parameters 
leads to strong resonant coupling and strongly perturbed 
vibrational wavefunctions where each diabatic compo- 
nent shows peaks at the four classical turning points of 
both potentials. As illustrated in the middle panel of 
Figf5J such a situation is potentially favorable for vibra- 
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TABLE I: Parameters of the spin-orbit coupling functions, cf. Eq. (J9]). 
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FIG. 5: (color online) Initial and target vibrational wavefunc- 
tions of the electronic ground state (upper panel) and effect 
of the spin-orbit coupling on the vibrational wavefunctions of 
the electronically excited states (medium and bottom panel). 



tional stabilization: The outermost peak of the vibra- 
tional wavef unction with a binding energy of 112 cm -1 
(corresponding to an absolute energy of 12625 cm- 1 ) 
leads to good Franck-Condon overlap with the initial 
state. The singlet component of the wavefunction shows 
a second peak at the outer turning of the 6 3 II state, 
R » 19 ao. This second peak will lead to better overlap 
with more deeply bound levels in the electronic ground 
state and could thus cause a speed-up of the stabilization 
dynamics toward shorter distances or less required pulse 
energy. Note that the Franck-Condon overlap reflects 
only the singlet component of the vibrational wavefunc- 
tions. However, due to the time-dependence of the stabi- 
lization process, the triplet component may play a role as 
well. If there is a dynamical interplay of pulse and spin- 
orbit coupling, the transfer efficiency can be much larger 
than predicted by static Franck-Condon overlaps p5j . 
Such a situation occurs for strong pulses and pulse du- 
rations comparable to or longer than the period of the 
singlet-triplet oscillations caused by the spin-orbit inter- 
action. Both conditions will be met by the optimized 
pulses presented below. 

In case 2, we have only modified the width of the dip in 
the off-diagonal spin-orbit coupling, P™. As in case 1, 
strong perturbations in the vibrational wavefunctions are 
observed, cf. the peaks at the two outer turning points in 
the singlet component in the lower panel of Fig. [SJ How- 
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FIG. 6: Optimal pulses (upper panel) and their spectra (lower 
panel) for the spin-orbit coupling cases 1 (left) and 2 (right). 



ever, the spin-orbit coupling is now expected to have a 
detrimental effect on vibrational stabilization where the 
wave packet shall be transferred from large to short dis- 
tances: Once the wavepacket comes close to the outer 
turning point of the upper adiabatic potential, R rs 21 ao, 
the resonant spin-orbit coupling will move part of the 
probability amplitude all the way out to the outer turn- 
ing point of the lower adiabatic potential, R rs 32 ao. 
Therefore in case 2, the spin-orbit coupling will poten- 
tially counteract the vibrational stabilization. 

The same three-step optimization procedure as in Sec- 
tion lllll has been followed: (i) optimization for T — 16 ps, 
(ii) compression in time to T — 4ps and subsequent 
re-optimization, (iii) determination of the minimal in- 
tegrated pulse energy with which a population trans- 
fer of better than 99% can be achieved. However, the 
guess pulses for step (i) were chosen such as to take 
the modified Franck-Condon factors into account. Fig- 
ure IH1 compares the optimal pulses and their spectra for 
spin-orbit coupling cases 1 and 2. Compared to Fig. [3] 
where the spin-orbit coupling in the excited state was 
completely neglected, the optimal spectra of Fig. [5] are 
broader, with additional spectral amplitude at small and 
large frequency components. However, as explained in 
Section Mil further calculations employing spectral fil- 
tering are necessary to determine whether these spectral 
features are artifacts of the optimization algorithm or 
whether they represent a true physical requirement that 
the optimal pulse has to fulfill. As seen in Fig. [6l the 
minimum optimization time to yield a population trans- 
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pulse energy than the estimate of 61 /j,J obtained with the 
two-state model of Section HID This is most likely due 
to the much larger state space that is explored by the 
optimization. 

The dynamics under the optimal pulses are analyzed 
in Fig. [7] by projecting the time-dependent wave packet 
onto the triplet and singlet components of the vibrational 
eigenstates. Overall the dynamics are very similar for the 
two spin-orbit coupling cases. A difference would be ex- 
pected mainly at the beginning of the pulse where the 
initial state is excited into levels of about 12630 cm -1 
energy of the electronically excited state. Inspection of 
the Franck-Condon map displayed in Fig. [5] (inset) shows 
that the resonant coupling leads to additional peaks com- 
pared to the model without spin-orbit interaction. These 
features are caused by the additional peak in the sin- 
glet component of the vibrational eigenfunctions near the 
outer turning point of the triplet potential, cf. Fig. [5] 
However, the projections of the wave packet within the 
first 1 ps do not reveal any substantial differences between 
coupling case 1 and 2. We therefore conclude that the 
modifications of the Franck-Condon map due to the spin- 
orbit interaction are not significant enough to influence 
the optimized stabilization dynamics, no matter whether 
the type of coupling is potentially favorable or poten- 
tially detrimental. The complete stabilization dynam- 
ics is rationalized in terms of the Franck-Condon map 
analogously to Section Hill i.e. it is determined by the 
main ridges of the Franck-Condon map. The only differ- 
ence between spin-orbit coupling cases 1 and 2 that can 
clearly be identified is the spread of population over the 
vibrational levels which is larger for case 2. 

To summarize the investigation of the influence of the 
spin-orbit interaction in the electronically excited state, 
the minimum optimization time is not affected while the 
required pulse energy is significantly increased compared 
to the model without spin-orbit interaction. Details of 
the spin-orbit interaction have only a minor effect on the 
required pulse energy and stabilization dynamics. 



fer of better than 99% is not affected by the spin-orbit 
coupling in the excited state. This is in accordance with 
the rationalization in terms of the time scale of the vi- 
brational dynamics on the electronic ground state, i.e. 
in terms of the timescales related to resolving the initial 
state and the target state, cf. Section UITl The mininum 
required pulse energy amounts to 140 /xJ for spin-orbit 
coupling case 1 and 180 /xJ for case 2. While the opti- 
mization target, i.e. population transfer of better than 
99%, can be achieved in both cases, case 2 which had been 
identified as potentially bad for the stabilization, requires 
more pulse energy. Both cases require substantially more 



V. OPTIMIZING TRANSFER FROM A 

SINGLET-TRIPLET SUPERPOSITION TO THE 

SINGLET GROUND STATE 



In heavy heteronuclear alkali dimer molecules it is pos- 
sible to transfer a vibrationally excited state that is in the 
lowest triplet state or in a superposition of the lowest 
triplet and singlet electronic ground state, to the rovi- 
bronic ground state [8]. In order to study this as an 
optimization problem, a four-state model of the KRb 
molecule is considered, 



H 4s (i) 



ff + V x ij,+ (R) Ae*(*) 

f + V a3s+ (R) 

A«;(*) o f + tWR) 

V jXe%{t) Wfg(R) t 



r 



o 

V b s n (R)-W™(R)J 



(10) 



that allows for transfer of molecules from the lowest 
triplet state to the singlet ground state due to the spin- 
orbit interaction in the electronically excited state. The 
lowest triplet state potential is taken from Ref. [37|, the 
other potential curves, dipole moments and spin-orbit 
coupling functions are constructed as described in the 
previous sections. In particular, the spin-orbit coupling 
cases 1 and 2 introduced in Section IIVI are employed. 
The initial state is taken to be a superposition of the vi- 
brational eigenfunctions of the a 3 E + lowest triplet state 
and X 1 S + singlet electronic ground state with 4.5 cm -1 
binding energy. The triplet (singlet) component carries 
70% (30%) of the population. The target state remains 
unchanged compared to the previous sections, i.e. the 
v = level of the X 1 S + singlet electronic ground state. 

Different laser polarizations need to be taken into ac- 
count - linearly polarized light, E n (t), for the singlet tran- 
sitions and circularly polarized light, £ CT (t), for the triplet 
transitions. This simply means that instead of Eq. © 
two equations for the two components of the field need 
to be considered where the dot products are evaluated 
for the corresponding components of the states. 

As explained above, a three-step optimization proce- 
dure is carried out in order to determine the minimum op- 
timization time and minimum pulse energies. Also for the 
four-state model, population transfer with an efficiency of 
better than 99% is achieved by the optimal pulses. While 
the optimization time remains unchanged at 4ps, the re- 
quired pulse energy is increased. It amounts to 270 fi3 
for 7r-polarization and 230 /iJ for er-polarization in spin- 
orbit coupling case 1, and to 300 /i J for 7r-polarization 
and 270 /iJ for er-polarization in case 2. As in Section HVl 
case 1 which is potentially favorable for the stabiliza- 
tion requires slightly less pulse energy than case 2. How- 
ever, the further increase of pulse energy compared to 
the three-state model of Section [IV] where pulse energies 
of 140 /iJ and 180 /xJ were obtained, is significant, in par- 
ticular in view of the fact that the molecule now couples 
to two polarization components. The high pulse ener- 
gies reflect the more difficult optimization problem that 
is considered here where the wave function needs to be 
changed qualitatively from a singlet-triplet superposition 
to a pure singlet state. 

The dynamics under the optimal pulses are analyzed 
in terms of the projections of the time-dependent wave 
packet onto the components of the eigenfunctions on the 
four electronic states in Fig. |H1 The dynamics of the sin- 
glet components is rather similar to those of the previous 
sections, cf. Figs. [3] and [7] Most of the triplet compo- 
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FIG. 9: (color online) Projection of the time-dependent wave 
packets onto the vibrational eigenstates of the electronic 
ground (bottom) and the coupled excited (top) states for spin- 
orbit coupling cases 1 (left) and 2 (right). 



nent of the initial state is converted to singlet compo- 
nents within the first half of the pulse, in particular in 
spin-orbit coupling case 1. In the potentially detrimen- 
tal case 2, more population resides in the triplet compo- 
nents than in case 1 and remains there throughout the 
pulse. The population dynamics should be compared to 
the Franck-Condon maps shown in Fig. 1101 The sin- 
glet dynamics start out in the upper right corner of the 
Franck-Condon maps. They follow the main ridges of the 
map but fan out as the dynamics roll down the ridge and 
levels with equally likele transitions to many levels are 
populated. The triplet dynamics follows the right-most 
ridge in the upper right corner of the Franck-Condon 
maps. It never jumps over to the left- most ridge; and 
in particular in case 1, the range between v' — 275 and 
v' = 50 of the 6 3 II state seems to be bridged via singlet 
dynamics. This is an indication for cooperative behavior 
between pulse and spin-orbit coupling in order to achieve 
the triplet-singlet transfer required by the optimization 
task: The dynamics travel down the right-most ridge of 
the Franck-Condon map due to Rabi cycling, i.e. due to 
interaction with the circularly polarized pulse. Then the 
spin-orbit coupling transfers most of the triplet popula- 
tion to the singlet channel, and a further decrease of the 
vibrational excitation happens in the singlet channels due 
to interaction with the Tr-polarized light. Such a coop- 
erative behavior can be expected since the optimization 
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time is much larger than the time scale corresponding to 
the spin-orbit interaction. In the potentially detrimental 
case 2, no indication for cooperative behavior between 
pulse and singlet-triplet oscillations is observed; and pop- 
ulation that was initially in the triplet channels remains 
there for a much longer time. 

To summarize this section, vibrational stabilization 
from a singlet-triplet superposition to a pure singlet state 
can be achieved with better than 99% efficiency. How- 
ever, even more pulse energy for both polarization com- 
ponents is required compared to stabilization of a pure 
singlet state. Depending on the details of the spin-orbit 
coupling function, cooperative behavior between pulse 
and spin-orbit coupling may or may not be observed. 



VI. CONCLUSIONS 

We have studied the vibrational transfer of KRb 
molecules from a level just below the dissociation limit 
to the vibrational ground state. Optimal control theory 
was employed to obtain the shaped laser pulses that drive 
this population transfer with an efficiency of 99% or bet- 
ter. As the main result, our calculations have yielded 
an estimate of the minimum time that is required for 
the vibrational transfer, i.e. the quantum speed limit for 
this process [46j, and an estimate on the required pulse 
energy. 

Our findings have confirmed that optimal control ap- 
proaches work as 'black-box' algorithms that provide so- 
lutions independent of the details of a given quantum ob- 
ject. Nevertheless, our results cannot straightforwardly 
be transferred to a laser pulse shaping experiment on 
cold molecules. Insufficient knowledge of the molecular 
structure and restriction of our model to a few electronic 



states prevent direct experimental application of the cal- 
culated pulses. This is a common phenomenon encoun- 
tered in the optimal control of complex quantum systems 
such as molecules [47| as opposed to atoms [48| or spin 
systems [49|. In principle, our model could be refined. 
For example, Ref. [12| discusses in detail how potential 
loss channels such as multiphoton ionization could be in- 
corporated; and more spectroscopy could be performed 
to obtained a better knowledge of the potential curves 
and non-adiabatic couplings. However, such refinement 
would miss the point of this study. Here, we come to 
a two-fold conclusion: On one hand, our study encour- 
ages optimal control experiments because solutions will 
be found no matter what are the specific details of the 
molecule. On the other hand, our study has clarified, 
as discussed below, the influence of the molecular struc- 
ture relevant for the vibrational transfer by successively 
increasing the complexity of the model. While of less im- 
portance in optimal control experiments based on feed- 
back loops, these findings are important for vibrational 
transfer and vibrational cooling using cw lasers, incoher- 
ent broadband light or adiabatic passage. 

First, we have addressed the role of the long-range be- 
havior of the excited state potential on the vibrational 
stabiliziation by comparing the KRb and Na2 molecules. 
From a time-dependent perspective one might expect 
the 1/R 3 potential of the homonuclear sodium dimer to 
be more favorable for the vibrational transfer since the 
larger slope of the potential speeds up the motion to- 
ward shorter internuclear distances. However, from a 
time-independent perspective, one might argue that the 
1/R 6 excited state potential of heteronuclear molecules 
yields better Franck-Condon overlap with the electronic 
ground state that also shows a 1/R 6 dependence at large 
internuclear distances. Comparing a two-state model, i.e. 
a model comprising of the singlet electronic ground state 
and a single excited state, for Na2 and KRb, we found 
that significantly less pulse energy is required for KRb. 
We therefore conclude that a 1/R 6 excited state poten- 
tial is more favorable for vibrational transfer than a 1/R 3 
potential, i.e. the spectroscopic perspective prevails over 
the dynamical one. 

The comparison of the KRb and Na2 molecules has 
also allowed us to identify what determines the minimum 
optimization time, i.e. the quantum speed limit for vibra- 
tional stabilization. The longest timescale in the problem 
that needs to be resolved by the optimal pulse is the vi- 
brational motion of the initial state. Taking comparable 
binding energies of the initial state for KRb and Na2, the 
difference in the vibrational periods is due to the mass 
of the molecules, and the minimum optimization time is 
smaller for the lighter molecule. 

Second, we have increased the complexity of the model 
for the molecular structure by taking spin-orbit coupling 
in the electronically excited state into account. The re- 
sulting model consists of three electronic states where the 
two electronically excited states exhibit a non-adiabatic 
coupling. In heavy alkali dimer molecules, the spin- 
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orbit interaction does not only modify the potentials at 
large internuclear separation, it may also cause a mix- 
ing of vibrational ladders affecting the complete vibra- 
tional spectrum. This effect has been termed 'resonant 
coupling' [2lM23| . Since the independence of the spin- 
orbit coupling function for KRb is not precisely known, 
we have adapted a parametrization developed by Berge- 
man et al. for RbCs [44J. In order to see which effect 
resonant spin-orbit coupling may have on the stabiliza- 
tion dynamics, we have modified the parameters of the 
coupling yielding a potentially favorable and a poten- 
tially detrimental case. While the resulting parametriza- 
tion may be far from the true spin-orbit coupling found 
in the KRb molecule, this approach allows us to iden- 
tify the maximum influence that the spin-orbit coupling 
may have on the stabilization dynamics. To our surprise, 
we found that in both the potentially favorable and the 
potentially detrimental coupling case, the pulse energy 
required to achieve population transfer to better than 
99% is significantly increased compared to the two-state 
model neglecting the spin-orbit interaction. This means 
that the increased size of the state space has a much 
larger effect on the optimization than a modification of 
the Franck-Condon factors underlying the dynamics. If 
compared amongst each other, the potentially detrimen- 
tal spin-orbit coupling case requires more pulse energy 
than the potentially favorable one. Overall, however, the 
details of the spin-orbit interaction seem to have only a 
minor effect on the stabilization dynamics. 

Third, the spin-orbit interaction in the electronically 
excited state allows for population transfer from a singlet- 
triplet superposition to a pure singlet level. In order to 
investigate this as an optimization problem, our mini- 
mal model consists of four electronic states, the singlet 
ground and lowest triplet state and the non-adiabatically 
coupled electronically excited states. Due to the symme- 
try of the electronic states, different polarization com- 
ponents of the laser field couple to transitions between 
the singlet and triplet channels. We found that the fur- 
ther increase in size of the state space as compared to the 
three-channel model results in even higher required pulse 
energies to achieve population transfer of 99% or better. 
For a shape of the spin-orbit coupling function that is 
potentially favorable to the vibrational transfer, we have 
observed indication for cooperative behavior between the 
pulse and the singlet-triplet transfer due to spin-orbit 
coupling. We have not seen any evidence for cooperative 
behavior in the case of potentially detrimental spin-orbit 
coupling. Correspondingly, the required pulse energy is 
larger for the potentially detrimental spin-orbit coupling 
case. Since the triplet-singlet transfer is explicitly part of 
the optimization problem in the four-state model, it is not 
surprising that details of the spin-orbit interaction play 
a somewhat larger role than for the three-state model for 
singlet-to-singlet vibrational population transfer. 

In summary, independently of the details of the molec- 
ular structure, we have found optimal pulses achieving vi- 
brational population transfer of KRb molecules to the vi- 



brational ground state with 99% efficiency or better. This 
highlights the power of the optimal control approach. 
However, as the complexity of the molecular structure 
is increased, the optimal laser fields need to carry more 
and more pulse energy. Each individual solution does of 
course depend on the details of the molecular structure, 
and we have analyzed the dynamics under the optimal 
pulse in terms of the underlying Franck-Condon maps. 

In the present paper, the initial state was taken to be a 
highly excited but pure state. Such a situation is encoun- 
tered for example if the molecules are created utilizing 
a Feshbach resonance [I3. Transfer to the vibrational 
ground state can then be achieved in a purely coherent 
process where the pulse absorbs the vibrational excita- 
tion energy of the molecule. The optimization task be- 
comes more involved if the initial state corresponds to an 
incoherent ensemble of vibrationally excited molecules. 
Such a situation occurs if the molecules are created by 
photoassociation followed by spontaneous emission [3j . A 
true cooling scheme is then required where the molecules 
can dispose of energy and entropy [l(| . While the present 
study was confined to vibrational stabilization, it may 
nevertheless shed light on the prospects for vibrational 
cooling. The Franck-Condon map shows a characteristic 
parabola whose distribution of weights and tilt determine 
whether the probability for heating is larger than that for 
cooling or vice versa. Similarly to the case of LiCs dis- 
cussed in Ref. [5fJ, the Franck-Condon maps presented 
above for Na2 and KRb reveal that vibrational cooling 
by optical pumping with a spectrally cut femtosecond 
laser pulse [1(| will not be successful. In order to prefer- 
entially cool instead of heat despite the Franck-Condon 
map, a more sophisticated approach than optical pump- 
ing would be required. Here again, optimal control can 
serve as the tool of choice. For a toy molecular model, 
optimally shaped laser pulses together with spontaneous 
emission have been predicted to yield a successful cooling 
scheme [5 If . 

Finally, we point out that we do not obtain an adi- 
abatic passage-like solution since this is not accessible 
within our current optimization approach [52l j . In fact, 
Stimulated Raman Adiabatic passage and related solu- 
tions formally require the infinite time limit [53|; and a 
finite optimization time needs to be fixed when numeri- 
cally solving the optimization equations. So on one hand 
optimal control is an extremely convenient tool that al- 
lows for solving very complex optimization problems. On 
the other hand, however, it is not always straightforward 
to translate physical considerations such as allowing for 
the adiabatic limit into mathematical prescriptions for 
the algorithm. Our future work is therefore dedicated to 
developing more versatile optimization algorithms. 
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